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Abstract 

Bell-type inequalities and violations thereof reveal the fundamental differences between standard prob- 
ability theory and its quantum counterpart. In the course of previous investigations ultimate bounds on 
quantum mechanical violations have been found. For example, Tsirelson's bound constitutes a global upper 
limit for quantum violations of the Clauser-Horne-Shimony-Holt (CHSH) and the Clauser-Home (CH) in- 
equalities. Here we investigate a method for calculating the precise quantum bounds on arbitrary Bell-type 
inequalities by solving the eigenvalue problem for the operator associated with each Bell-type inequality. 
Thereby, we use the min-max principle to calculate the norm of these self-adjoint operators from the maxi- 
mal eigenvalue yielding the upper bound for a particular set of measurement parameters. The eigenvectors 
corresponding to the maximal eigenvalues provide the quantum state for which a Bell-type inequality is 
maximally violated. 

PACS numbers: 03.67.-a,03.65.Ta 

Keywords: Bell-type inequalities, quantum probabilities 



*Electronic address: jsfilipp@ati.ac.at| 



^Electronic address: lsvozil@tuwien.ac.atl URL: |http : //tph . tuwien . ac . at / ~ svozil| 



1 



I. INTRODUCTION 



One of the most puzzling features of quantum mechanics is the violation of so-called Bell-type 
inequalities representing a cornerstone of our present understanding of quantum probability the- 
ory 111]. As pointed out by John Bell |Q] such a violation, as predicted by quantum mechanics, 
requires a radical reconsideration of basic physical principles like the assumption of local realism. 
However, Bell-type inequalities have already a long tradition dating back to George Boole's work 
on "conditions of possible experience" i^m, dealing with the question of necessary and sufficient 
conditions on probabilities of logically interconnected events ll42ll . Take for example the state- 
ments: "The probability of rain in Vaxjo is about 80%" and "The probability of rain in Vienna is 
90%". Nobody would believe that the joint probability of rain in both places could be just 10% 
— the claim that the joint probability is very much lower than the single probabilities is appar- 
ently counterintuitive. The question remains: Which numbers could be considered reasonable and 
consistent? Boole's requirements on the (joint) probabilities are expressed by certain equations or 
inequalities relating those (jointVprobabilities. 

Since Bell's investigations ||2|, |5|] into bounds on classical probabilities and their relation to 
quantum mechanical predictions, similar inequalities for particularphysical setups have been dis- 



cussed in great number and detail (see for example Refs. BaLZllSlla])- Furthermore, violations of 
Bell-type inequalities, as pred icted b y q uantum mechanics, have been experimentally verified in 
different areas of physics lioi nl[l2[ |l3l to a very good degree of accuracy. 

However, whereas these bounds are interesting for an inspection of the violations of classi- 
cal probabilities by quantum probabilities, the issue of the validity of quantum probabilities and 



their experimental verification is completely different. Recently, Bovino et al. 11141] conducted an 
experiment based on numerical studies by the current authors lisl and triggered by a proposal 
of Cabello 1 3] to verify bounds on quantum probabilities depending on a particular choice of 
measurements. 

In what follows we shall present analytical as well as numerical studies on such quantum 
bounds allowing for further experimental tests of different kinds of Bell-type inequalities. 
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II. CORRELATION POLYTOPES 



At first we shall start from a geometrical derivation of bounds on classical probabilities given 
by linear inequalities in terms of correlation poly topes 1 17, l^ll^lzQ. 21\- Considering an ar- 
bitrary number of classical events ai,a2, ■■ .,an one can assign to each event a certain probability 

jci, p2, ■ ■ ■, Pn and probabilities pi2, ... for the joint events ai Ciaj, These probability values 

can be collected to form the the components of a vector p = {pi,P2t ■ ■ , Pm Pii^ - ■ where each 
Pi, pij {i,j= 1, . . . ,n) can take values in the interval [0, 1]. Since the events ai, a2, ... are as- 
sumed to be independent, each single probability pi can also take its extremal value or 1 and 
the vectors comprising all possible combinations of extremal values (p, = 0, 1 and pij = PiPj) can 
be regarded as rows of a truth table; with the symbols "0" and "1" corresponding to "false" and 
"true," respectively. 

Any classical probability distribution; i. e., any vector p, can be represented as a convex sum 
over the extremal probability distributions given by the row entries of the truth table. It can there- 
fore be regarded as some point p EC where C = conv{K) is a convex poly tope defined by the set 
of all points that can be written as a convex sum extending over all vectors associated with row 
entries in the truth table. More formally, 



conv{K) 

with 





(1) 



{xi,X2,...,X2«} = { {ti,t2,...,t„,tjy,...) I tie{0,l},i=l,...,n}. (2) 

Here, the terms txty, . . . stand for arbitrary products associated with the joint propositions which 
are considered. Exactly what terms are considered here depends on the particular physical config- 
uration. 

In a next step towards the linear inequalities sought one utilizes the Minkowsky-Weyl repre- 
sentation theorem p. 29] stating that every convex polytope in Euclidean real space has a dual 
description: either as the convex hull of its extreme points - in our case the rows of the truth table 
- or as the intersection of a finite number of half-spaces. Each half space can be described by a 
linear inequality. To obtain the inequalities from the vertices one has to solve the so-called hull 
problem. These inequalities coincide with Boole's "conditions of possible experience"; i. e., they 
constitute the bounds of classical probabilities. The set of inequalities obtained is maximal and 
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complete, as no other system of inequalities exists which characterizes the correlation polytope 
completely and exhaustively. 

For particular physical setups these inequalities correspond to Bell-type inequalities. Therefore 
correlation polytopes provide a constructive way of finding the entire set of Bell-type inequalities 



for a given physical configuration ||2 



2411 . although from a computational complexity point of 



view the problem remains intractable ll2f 

As an example, we consider the derivation of the well known Clauser-Horne inequality 
Given a source emitting pairs of correlated spin- 1/2 particles either in the positive or in the neg- 
ative 3;-direction, the spin of both particles can be measured in arbitrary directions perpendicular 
to the propagation direction; i. e., restricted to the x-z plane. Implementing two measurement 
directions on each side labeled by the angles a, P for the particle propagating in the negative y 
direction (left hand side) and y, 5 for the particle propagating along the positive y axis (right hand 
side), we obtain the probabilities for measuring "spin-up" for each particle and measurement di- 
rection PaiPpiPyiPs- The joint probabilities for measuring "spin-up" on the left while measuring 
"spin-up" on the right in coincidence - but in general with different measurement directions - 
are denoted by Pay^Padi ■ ■ ■ The probability distribution vector for this situation is consequently 
P = (Pa: P^: Pj: P8: Pay: P^y: P^8: Pad) ^^d thc truth table (comprising the extremal probabilities) 
consists of 2^ = 16 rows by inserting pa^ i?p, i?y, p§ G {0, 1}. The corresponding polytope is eight 
dimensional. By solving the hull problem, which for this simple setup can easily be done, we ob- 
tain inequalities like < pa, Pay < ^, Pa + Py — Pay < 1 ; and in the similar manner for p^, py, p^. 
Additionaly, inequalities of the form 

-'^ < Pay + Pa8+Ppy-P^?,-Pa-P8<0 (3) 

also represent bounds of this correlation polytope. The inequality ©, termed Clauser-Home 
(CH) inequality, and the inequalities containing all permutations of the parameters, are violated 
by quantum theory for particular choices of the angles and for specific quantum states. They 
constitute therefore a demarcation criteria between classical and non-classical probabilities, such 
as the ones encountered in quantum theory. 
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III. VIOLATION OF BELL-TYPE INEQUALITIES 



Similar to the bounds on classical probabilities given by the Bell-type inequalities, there exist 
bounds on quantum probabilities which will be the subject of the following discussion. There have 



been investigations in the analytic aspects of bounds on quantum probabilities, most prominently 
by Tsirelson 112 IL 12611 and recently by others in Refs. but also numerical 1 15] and 



experimental lll4ll test have been performed. The quantum probabilities do not violate Bell-type 
inequalities maximally 1291 BoHbiI . Take, for example, the well known Clauser-Horne-Shimony- 



ity B 



Hok (CHSH) inequality 

|£(a,Y)-f£(|3,Y)+£(|3,5)-£(a,5)|<2, (4) 

where £'(/j,v) denotes the correlation function for two particle correlations with possible values 
in the interval [—1,1] when measuring their spin/polarization in coincidence along the directions 
/J and V, respectively. The global limit for a quantum violation of this inequality is 2\f2. I^llls^; 
quantum theory does not allow a higher value, no matter which state and which measurement 
directions are chosen. However, in principle, the four terms on the left hand side of Eq. Q could 
be set such that a value of 4 can be obtained by appropriate choices of ±1 for the correlation 
functions. 

Popescu and Rohrlich investigated the case where "physical locality" is assumed without 
referring to a specific physical model (such as quantum mechanics), whether realistic or not. In 
this context, "physical locality" means that the marginal probabilities for measuring an observable 
on one side should be independent of the observable measured on the other side, which is a natural 
assumption for a Lorentz invariant theory. The maximal value of the left hand side of Eq. © 
has been shown to be 4 as well, which is beyond the quantum bound 2\f7. and we can conclude 
that quantum mechanics does not exploit the whole range of violations possible in a theory con- 
forming to relativistic causality. Still, in our opinion, the nagging question remains why quantum 
mechanics does not violate the inequality to a higher degree. 

In what follows, we will restrict our attention to the simpler task to explore the quantum bounds 
on violations of Bell-type inqualities for particular given measurement directions and arbitrary 
states. It turns out that the equations for the analytic description of the quantum bounds can be 
derived by solving an eigenvalue problem. Intuitively it cannot be expected that it is feasible to 
achieve a maximal violation of some inequality for any set of measurements just by choosing a 
single appropriate state. 



The quantum mechanical description of the physical scenario discussed above involves spin 
measurements represented by projection operators 

F(0) = 1(12 + 5(0)), (5) 

(COS0 sin0 \ 
, denoting the direction of measurement in the x-z plane, and 
sin0 — COS0 / 

I2 standing for the two-dimensional identity matrix. For an even more general description we 
would have to take all possible two-dimensional projection operators into account, corresponding 
to measurements in arbitrary directions. As this generalization is straightforward and does not 
lead to any more insight, we will work with this restricted set of measurements parameterized in 
Eq. ©. 

F(0) acts on one of the two particles. This implies that we have to choose a tensor product 
of two Hilbert spaces to represent the state vectors corresponding to possible state configurations; 
i. e., 9{ = H2- The representation of a single-particle measurement in 0-C is then 

^(e) = ^(l2 + o(0))®l2 or l2®i(l2 + o(0)) (6) 

for a measurement on the particle emitted in the negative y-direction {0-C\), or in the positive y- 
direction respectively. Two-particle measurements are implemented by applying i^(0) on 
both !}{\ and iHr, i. e., 

^(0,0O = ^(l2 + o(0))®^(l2 + 6(0')), (7) 

corresponding to a measurement of the joint probabilities. This setup can easily be enlarged to sys- 
tems comprising more than two particles by the tensor product of the appropriate Hilbert spaces, 
but for the sake of simplicity we will restrict ourselves to bipartite systems. 

The general method for obtaining the quantum violations of Bell-type inequalities is then to 
replace the classical probabilities by projection operators in Eqs. ( I6I7I) in a certain Bell-type in- 
equality to obtain the Bell-operator, which is a sum of projection operators. In the case of the CH 
inequality one obtains 

-1 < ^(ea,0y) +^(0a,05) +^(e(3,0Y) -q{%%)-q{Qa) "^(65) < 0. (8) 
In a second step one calculates the quantum mechanical expectation values by 



(^(0))=Tr[W^(0)], 
6 



(9) 



where V7 is a positive definite, Hermitian and normalized density operator denoting the state of 
the system. For some W and set of angles {0a, 6(3,67, 65} one obtains a violation of a classical 
inequality. 

In general the Bell-operators can be written in the form 

0= £ bi^i^,„i^Pi,®Pi^®...P,^, (10) 

with real valued coefficients bi^i2,„ij^. Here A^^ is the number of particles involved and the Pj axe 
either projection operators denoting a measurement on particle i or the identity when no mea- 
surement is performed on the z'-th particle. Since (A 5)^ = = A0B and {A + BY = 
A^ +B^ = (A + B) for arbitrary selfadjoint operators A,B, the Bell-operator O is also self-adjoint 
with real eigenvalues. However, the eigenvalues of O cannot be deduced from the eigenvalues of 
the constituents Pi^ ® /'ij ® . . - Pi,, in the sum in Eq. (flOb since these are not commuting in general 
and therefore are not diagonalizable simultaneously. 

One can make use of the min-max principle §90], stating that the bound of a self-adjoint 
operator is equal to the maximum of the absolute values of its eigenvalues. Thus, the problem of 
finding the maximal violation possible for a particular choice of measurements can be solved via 
an eigenvalue problem. The maximal eigenvalue corresponds to the maximal violation and the 
associated eigenstates are the multi-partite states which yield a maximum violation of the classical 
bounds under the given experimental (parameter) setup l^^. 

For a demonstration of the method let us start with the trivial setup of two particles measured 
along a single (but not necessarily identical) direction on either side. The vertices are (pi , pi.pii ~ 
PiPi) for Pi-,P2 G {0, 1} and thus (0,0,0), (0, 1,0), (1,0,0), (1, 1, 1 ); the corresponding face (Bell- 
type) inequalities of the polytope spanned by the four vertices are given by pi2 <P2,Q< Pn < 1> 

andpi-fjc»2-pi2 < 1- 

The classical probabilities have to be substituted by the quantum ones; i.e.. 

Pi ^ ^i(e) = i[l2 + o(0)]®l2, 

P2 - ^2(e) = l2®^[l2 + a(0)], (11) 

pn ^ ^i2(0,0O = 5[l2 + o(0)]®i[l2 + o(0O]. 
It follows that the self-adjoint transformation corresponding to the classical Bell-type inequality 
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{pi + Pi-pii < 1) is given by 



On (O,0)=^i(O)+<5r2(e)- ^12(0,0) 



/l 
1 








\ 




cos^l ^ 



si§e sin2 § 



\0 



(12) 



The eigenvalues of On are and 1, irrespective of 0, the maximal value of On predicted by 
the min-max principle does not exceed the classical bound 1. 

Now we shall enumerate analytical quantum bounds for the more interesting cases comprising 
two or three distinct measurement directions on either side yielding the quantum e quiv alents of 
the Clauser-Horne (CH) inequality, as well as of the inequalities discussed in IibI 13511361 . 

For two measurement directions per side, we obtain the operator O22 based on the CH- 
inequality [Eq. Q] upon substitution of the classical probabilities by projection operators: 



O22 (a, P, Y, 5) = (a, y) + <5ri4 (a, 5) + ^23 (P, y) - g24(P, 5) - ^1 (a) - ^3 (y) ■ 
The eigenvalues of the self-adjoint transformation in ([T3t are 

Xi,2,3,4(a,|3,Y,5) = ^ ( ± Vl ± sin(a - (3) sin(Y- 5) - l) , 



(13) 



yielding the maximal bound IIO22II = max,=i.2,3,4?t/. The eigenstates corresponding to maximal 
violating eigenstates are maximally entangled for general measurement angles lying in the x-z- 



plane 



The numerical simulation of the bounds of the CH-inequality is based on the generation of 
arbitrary bipartite density matrices W; i. e., 4 x 4 Hermitian positive matrices with trace equal to 
one. Since one can write a Hermitian positive matrix W as the square of a self-adjoint matrix, W ~ 
B^. The normalized matrix W' = VF/Tr[W] can thus be explicitly parameterized by 16 parameters 
bi,b2,...,bie e R; i. e., 

/ bi bs+ibe ^11 +i^i2 ^15 +i^i6 \ 
bs-ibe bi bq+ib^ Z?i3+i^i4 
bii—ibi2 bj—ibs b^ bg + ibio 
\bi5-ibi(, bi^-ibu bg-ibio b^ J 

For a particular choice of projection operators, one can then generate random states W' in order 
to find the maximal violation possible for the current set of projection operators. In Figure [H 



1 



i:=ibi+2i};L,bj 



(14) 
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both the analytic and the numerical bounds are depicted for measurement directions a = 0, [3 = 
20, y = and 5 = 30 dependent on a single parameter G [0,7r]. In addition, the well-known 
maximal violation for the singlet-state at = 7r/4 and = 37i:/4 is drawn. 




FIG. 1 : Numerical simulation of the bounds of the CHSH-inequality 



The extension to three measurement operators for each 



non-equivalent inequality (with respect to symmetries) ||3 



3, 



particle merely yields one additional 



3^ 



h-i = PU + Pl5 + Pl6 + P24 + P25 - P26 + P34 - P35 - Pi -2p4-p5 <0 (15) 

among the 684 inequalities representing the faces of the associated classical correlation poly- 
tope. The associated operator for symmetric measurement directions is given by 

O33 (0, 0, 20, 0, 0, 20) = ^14(0, 0) + ^15 (0, 0) + ^16(0, 20) + q24iQ, 0) + <?25 (6, 0) - ^26(0, 20) + 

+^34(20, 0) - ^35(20, 0) - ^i(O) - 2<5r4(0) - <?5(e) 

/-4sin20 \ 

-5-2cos0-3cos20 + 2cos30 4cos2| 2sin0 + 3sin20-2sin30 

4cos2| -2(3 + cos20) -2sin0 

\ 2sin0 + 3sin20-2sin30 -2sin0 2sin2 f cos^ |(4cos0 - 3) / 

(16) 

in the Bell basis = {(1,0, 0,0)^(0, 1,0, 0)^(0, 0,1, 0)^(0, 0,0,1)^} 

with \\\f^) = 1/^2(101) ± |10)) and = 1/V2(|00) ± |11)). In this basis, O33 can be decom- 
posed into a direct sum of a one-dimensional and a three-dimensional matrix O33 = o j ©03, thus 
simplifying the calculations of the real eigenvalues. By using the Cardano method 11371] . these can 
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be calculated to be 



^3,4 



-sin^e 

-2 



M cos , 

' 3 3 

cos - ± sin - 
3 3J 



b 
3' 



(17) 



Here, m = l/9(3c — b^) and cos^ = ^{9bc — 2Z?^ — lid)/ (ma/H) where = — Tro3, c = 
l/2^Tr^03 — Trojj, J = — det03. (For convenience we have omitted the dependencies on 0.) 
In Figure 121 the eigenvalues are plotted as functions of the parameter 0. The maximum violation 




FIG. 2: Eigenvalues of O33 in dependence of the relative angle 6. 
of 1/4 is obtained for = 7l/3 with the eigenvector corresponding to 

|^max) = ^r) + ^l¥+)- (18) 

l^max) is maximally entangled, but in contrast to the CH-inequality, this is in general not the 
case for eigenstates corresponding to the maximal eigenvalue at 7^ 7r/3. 



IV. RELATION TO EXPERIMENTS 

The analytical quantum bound of the CH-operator has been enumerated by Cabello as well 

^ n 

as by the current authors [28] and experimentally verified by Bovino et al. II14I1 using polarization- 
entangled photon pairs. The ansatz of Cabello for the experimental realization made use of the fact 
that the eigenstates leading to maximal violations are maximally entangled. Thus when applying 
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a unitary transformation U (0) of the form 

f/(e) = 



sinG — COS0 
COS0 sin0 



onto an initial state |\|/ ) = -^(|01) — |10)), one obtains all maximally violating states for different 
values |\^(0)) = f/(0) ® hIV")- 

However, in the case of O33, this scheme has to be extended, since the maximal violating states 
are not maximally entangled in general. Such states cannot be created from maximal entangled 
initial states by a local unitary operation t/2x2 ^ SU{2) iS>SU{2), since such a factorized trans- 
formation does not change the degree of entanglement. To obtain states constituting the quantum 
bounds, one has to apply unitary transformations U4 E SU (4) to the initial state comprising also 
non-local operations which cannot be written as a tensor product of two unitary single-particle 
operators. 

A simplification for an experimental verification of the quantum bounds of Bell-type inequali- 
ties is due to the fact that maximal violating states are pure. Therefore, it is sufficient to generate 
initial states with variable degree of entanglement. Utilizing the Schmidt-decomposition, which is 
always possible for a bipartite state, one can write any pure state in the form |\|/) =Y.k'^k\ki) \k2) 
where \kj) are orthonormal basis states for particle 1 and 2, respectively, and Efc^i = 1- The 
weights of the X^s, are a measure of the degree of entanglement comprising the special cases 
where = ?i2 = 1/ V^, for a maximally entangled state and ?ii = 0, = 1 (or vice versa) for 
a separable state. Having a source producing such states in a particular basis one can obtain all 
other pure states by applying a local unitary operation 1/2x2 ^ SU(2) x SU{2). Appropriate pho- 



ton sources have been suggested for example by White et al. 11381] and Barbieri et al. [39] and 
could therefore be used to trace the bounds on arbitrary bipartite Bell-type inequalities in the same 
manner as in the experiment of Bovino et al. [14]. 



V. CONCLUSION 



In conclusion we have shown how to obtain analytically the quantum bounds on Bell-type in- 
equalities for a particular choice of measurement operators. We have also presented a numerical 
simulation for obtaining these bounds for the CH-inequality. We have provided a quantitative 
analysis and derived the exact quantum bounds for bipartite inequalities involving two or three 
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measurements per site. The generalization to an arbitrary number of measurement parameters is 
straightforward as the dimensionality of the eigenvalue problem remains constant. For more than 
two particles the dimension of the matrix associated with a Bell-type operator increases exponen- 
tially. However, one may conjecture that such matrices can be decomposed into a direct sum of 
lower dimensional matrices. 

In the context of this conference we also believe that the analytic expressions of the quan- 
tum bounds could serve as consistency criteria of mathematical models proposed to show that 
a violation of Bell-type inequalities does not necessarily imply the absence of a possible local- 
realistic theory from the logical point of view. It is claimed that violations can be achieved without 
abandoning a local and realistic position assuming for example time-dependencies of the random 
parameters |40], or "chameleon" effects ll4li . Still, any appropriate model has to be in accor- 
dance with quantum mechanics not only qualitatively, but also quantitatively, and hence should 
reproduce also the "fine structure" of the quantum bounds as discussed above. 

Finally, although there is no theoretical evidence for a stronger-than-quantum violation what- 
soever, its mere possibility justifies the sampling of the fine structure of the quantum bounds from 
the experimental as well as the theoretical point of view in order to understand and verify the 
restriction imposed by quantum theory. 
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